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We performed light optical diffraction experiments on a nanoparticle-polymer volume holographic 
grating in an angular range including also far off-Bragg replay A comparison of three diffraction 
theories - on the same level of complexity - with our experimental results shows that the dynamical 
theory of diffraction and the first-order two-wave coupling theory using the beta- value method fit 
the data very well. In contrast, the prevalent two- wave coupling theory using the K- vector closure 
method yields a poor fit with an order of magnitude worse mean squared error. These findings must 
be considered for accurate determination of coupling strength and grating thickness. 



I. INTRODUCTION 

It seems that after about eighty years on theories of 
diffraction from periodic structures everything should 
have been said and done. This is even more true since 
Moharam and Gaylord in a series of papers published 
a rigorous theory of diffraction (rigorous coupled wave 
analysis, RCWA), which covers nearly the entire scope of 
cases that have ever been relevant [TUTU]. However, for 
practical purposes approximate theories are still around, 
which are used instead. While the optics community has 
strongly opted for the coupled-wave approach of Kogel- 
nik (K- vector closure method, KVCM) [11] and much less 
for the one introduced by Uchida (Beta- value method, 
BVM) [12] , the neutron and X-ray communities are using 
solely the dynamical theory of diffraction (DDT) [T3HT9] 
or Darwin's early variant [20] [21] . Despite the fact that 
the discussion on diffraction theories for volume gratings 
might sound outworn, we will show that it is possible to 
experimentally single out a set of analytic theories which 
is significantly superior to others. 

We identified three problems in the theories when 
treating the off-Bragg regime: The first is, that the 
KVCM misses to correctly account for energy conserva- 
tion, the second is that KVCM and BVM only treat the 
half-space case, i.e., a single boundary. Thus boundary 
conditions for a parallel slab are not properly included 
and yield an ambiguity. Finally, with the DDT, which 
in principle yields the exact solution to the wave equa- 
tion, typically an approximation of the dispersion surface 
for small off-Bragg conditions is performed in literature 
[I3j [I9j [22] , which leads to hyperbolic dispersion surfaces 
and is called hyperbolic approximation in what follows. 
It is stated by Syms and Soly mar that the differences be- 
tween the amplitudes derived from each of the theories 
are small [23 . However, for geometrically thin gratings 
with large coupling coefficients as for our samples devia- 
tions in the far off-Bragg regime are observable. 



martin. fally@univie.ac. at 



The result of this work is that the first order two- wave 
coupling theory employing the BVM as well as the DDT 
[22] fit the experimental data extremely well, provided 
that we refrain from using the hyperbolic approximation 
for the latter. The KVCM fits the experimental data 
much worse. 

The paper is organized as follows: we start with a 
derivation of the relevant equations for the two- wave cou- 
pling theories KVCM and BVM as well as for the modal 
DDT without the hyperbolic approximation. Then the 
experiment and the corresponding results are shown to- 
gether with fits to each of the theories. We also compare 
the theories using proper approximations, discuss the im- 
plications and end with a conclusion. 



II. DIFFRACTION FROM A SINUSOIDAL 
VOLUME GRATING 

In what follows we summarize the different approaches 
to solve the wave equation for sinusoidal transmission vol- 
ume gratings and give the results for the diffraction am- 
plitudes with a particular emphasis on the far off-Bragg 
regime, which is subject to our experiments. 

We start with the simplest set of assumptions. Space 
is divided into three distinct regions: the input region 
(free space), the grating region (periodic material), and 
the output region (free space). A sketch of the geome- 
try is provided in Fig. [I] The grating be lossless, one 
dimensional (modulated only along the x direction), si- 
nusoidal, isotropic, and of the pure phase-type. Thus it 
can be described by 

n{x) = no + ni cos (Kx), (1) 

where no is the mean refractive index of the material un- 
der investigation, n\ the refractive-index modulation and 
K the spatial frequency of the grating. The scalar wave 
equation (Helmholtz equation) , for the polarization state 
perpendicular to the plane of incidence (s or H mode po- 
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FIG. 1. Illustration of the geometry and the boundary condi- 
tion. ko,qRi qs are the wavevectors of the incident wave in the 
input region, the forward diffracted wave, and the diffracted 
wave in the grating region. K, N,<B,6b denote the grating 
vector, the sample surface normal unit vector, the angle of 
incidence (external) and the Bragg angle (in the medium). 



larization) in such a medium is 

V 2 + (n(x)ko) 2 } E(x,z) = 0, 



(2) 



where ko = 27r/A with A the free space wavelength of 
light and E(x, z) the normalized electric field amplitude. 
It was shown that an exact solution of the wave equa- 
tion can be obtained by solving an infinite number of 
coupled linear differential equations of first order under 
appropriate boundary conditions pQ (RCWA). In the rig- 
orous treatment of the problem the coupled wave analysis 
and the modal approach are completely equivalent, i.e., 
lead to identical results [2]. For thick volume holograms 
discussed here only two diffracted beams of considerable 
field amplitudes exist at the same time, e.g., the first 
and zeroth diffraction order [23] . In addition, the ampli- 
tudes strongly depend on how well the Bragg condition, 
2/3 sin Ob = K, is fulfilled. Here Ob is the Bragg angle 
and /3 = korio the propagation constant in the material. 
A commonly defined off-Bragg parameter [TT] is 



i? = if sin 



K 



(3) 



where is the angle of incidence. Our interest is to com- 
pare three common diffraction theories with experimen- 
tal data obtained for a wide range of angles, particularly 
in the far off-Bragg regime. In contrast, usually only 
the very vicinity of the Bragg peak is considered and a 
linearization of d is performed (see also appendix |B| , a 
simplification which is not applicable in our case. In what 
follows we give the relevant equations for each of the three 
theories that cover also the far off-Bragg regime. 



A. First-order two-wave coupling: Kogelnik's 
approach (KVCM) 

The basic idea of Kogelnik's coupled wave approach is 
to solve Eq. ([2| by the ansatz 

E(x, z) = R(z) exp (iqn • x) + S(z) exp (iqs ■ x), (4) 

i.e., the sum of two waves whose amplitudes vary when 
propagating through the sinusoidal refractive index pat- 
tern. By inserting Eq. Q into Eq. Q two coupled dif- 
ferential equations for the amplitudes R(z),S(z) result 
that can be solved for appropriate boundary conditions. 
In most cases, second order derivatives can be neglected 
because the amplitudes are slowly varying functions as 
compared to the exponentials in Eq. Q (slowly varying 
envelope approximation). This simplifies the system of 
differential equations to first order on the expense that 
the boundary conditions are to some extent ambiguous 
(we end up with a half-space case instead of a slab, i.e., 
the output region is not considered). Coupling of the 
waves arises via the Bragg condition, which relates the 
wavevector qs of the diffracted beam to the wavevec- 
tor qji of the forward diffracted - frequently called also 
'transmitted' [25 - beam. At this point Kogelnik ad hoc 
introduces the relation 



qs = q~k± K. 



(5) 



While this is a correct choice from a mathematical point 
of view [26] , it gives physically meaningful results only if 
the Bragg condition is fulfilled exactly. When going off- 
Bragg, Eq. ([5| predicts either a wavelength change (sort 
of 'inelastic' scattering) that does, in fact, not occur, or 
a refractive index that strongly depends on the off-Bragg 
parameter. Furthermore, the direction of the diffracted 
beam's wavevector is not correctly predicted for the off- 
Bragg case [27]. A wavevector diagram for the off-Bragg 
case using Kogelnik's choice (KVCM) is shown in Fig. [2j 
The diffraction efficiency for the first order, the quantity 
measured in our experiment, is given by 
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with n = riiir/X the coupling coefficient, d the hologram 
thickness, cr = cosO and * denotes the complex conju- 
gate. 



B. First-order two- wave coupling: Uchida's 
approach (BVM) 

The only difference between Uchida's and Kogelnik's 
approach is the choice of the diffracted wave vector as 



qs = qk±K + MoN. 



(9) 



3 




diffraction efficiency 



FIG. 2. Wavevector diagram for the KVCM (Kogelnik) in off- 
Bragg position with f3 = |<f#| / \qs\ according to Ref. [TJ. 
N the surface normal unit vector. 




FIG. 3. Wavevector diagram for the BVM (Uchida) in off- 
Bragg position with f3 — \qji\ = \qs\ according to Ref. [12] . 



The BVM ensures energy conservation, i.e., \qs\ = \c[r\ = 
/?, by introducing a phase-mismatch parameter Ago- 
The resulting wavevector diagram for the off-Bragg case 
(Ago 7^ 0) can be seen in Fig. [3j From geometrical rea- 
soning the phase mismatch amounts to Ago = P(cr-cs) 
with c^- 
ter X = 



= y^cos 2 6 + 2#//3 = cr^/1 + 2X and a parame- 
$/((3c 2 R ), which is convenient for approximations 
discussed later (Appendix [B| . Solving Eq. Q with the 
wavevector choice of Eq. (pi) yields for the first order 
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These equations that have been derived, e.g., in Ref. [28] 
look quite similar to the ones obtained for Kogelnik's 
approach. Uchida 's seminal paper already suggested the 
BVM [12] . However, it was somewhat hidden by the 
second important topic addressed: the attenuation of the 
grating modulation along the sample depth. The latter is 
not included in the present treatment in order to enable 
a direct comparison of the theories. 



C. Two- wave modal approach: Dynamical theory 
of diffraction (DDT) 

The modal approach to solve Eq. Q is well known 
in solid state physics for electrons in a periodic potential 
(band structure, e.g., [29, 30 J, and also for diffraction of 
X-rays (e.g., [El HS1 113 EH) or neutrons [13 by crystal 
lattices. In volume holography it has been much less 
prominent [22, 32, 33 . The reason might be that it can 
be usefully applied only for highly idealised gratings [22] . 
In the modal approach the solution is taken rigorously in 
the form of a sum of m Bloch waves (eigen-modes) with 
a grating periodic amplitude function 



kx 



(13) 
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Fourier transforming Maxwell's equations and inserting 
Eqs. ( p~3] ) and yields a system of coupled algebraic 
equations for the amplitudes S^ 1 . Details can be found 
in Appendix |A] and Ref. [22]. For sinusoidal volume holo- 
grams we know that only four of the amplitudes - two for 
each eigen-mode - have appreciable magnitude and thus 
a linear system of two equations for each eigen-mode re- 
mains. To obtain consistent and non-trivial solutions to 
the systems, conditions for the magnitudes of the per- 
mitted wavevectors in the grating region arise: 



(2/? K ) 2 = (/3 2 
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(14) 



The latter is the decisive equation in the DDT and repre- 
sents the so-called dispersion surface, i.e., the surface of 
permitted wavevectors. The dispersion surface and the 
wave vectors of the eigen-modes in the grating region are 
shown in Fig. [4] These permitted wave vectors can eas- 
ily be found geometrically by following the three steps 
below: 
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1. 



3. 



Starting from the origin of the reciprocal space 
draw the incident wavevector with length fco at the 
(external) angle of incidence G. 

Phase matching at the boundary requires that 
the tangential compontents of the wavevectors are 
identical, i.e., fco, a; = Qr,x- Thus draw a line along 
the sample surface normal N and find the points of 
intersection A, B with the dispersion surface (typi- 
cally called 'tie-points', two in the transmission ge- 
ometry) 

The permitted wavevectors of propagation are 
now constructed by forming the vectors Alt = 



a* K 



AO) 



4 1} 



q£',±?K =qs 

The circles with radius j3 around the reciprocal lattice 
points 0, K represent the surface of permitted wavevec- 
tors in the medium without a grating viz. k = 0. They 
approximate the permitted wavevectors of propagation 
in a grating very well unless in the very vicinity of the 
Bragg condition. This latter situation is shown in the 
closeup of the dispersion surface in the lower part of Fig. 

m 

The modulus of the permitted wave vectors of propa- 
gation in the grating region are obtained by using Bragg's 
condition q$ — ± K and solving Eq. (14) for q^ 



-+(m) 
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(15) 



This means, that for each propagation direction un- 
der consideration (q = 0,K) there are two eigen-modes 
m = 0, 1 with corresponding propagation constants viz. 
energies or - as initially pointed out by Ewald - a sort of 
birefringence. 

Usually, the region in the vicinity of the Bragg condi- 
tion is approximated, i.e., the spheres around the recip- 
rocal lattice points 0, K with radius j3 are approximated 
by planes (asymptotic gray lines in Fig. |4| and the dis- 
persion surfaces are hyperbolae. We disregard this sim- 
plification so that the theory can be applied for the far 
off-Bragg region, too. For neutrons the corresponding 
equations for the diffracted intensities were established 
in Ref. [34j [35] . The amplitudes are determined by 
taking into account the boundary conditions and can be 
expressed by the grating parameters. Then the first or- 
der diffraction efficiency according to the DDT takes the 
form 



where Aq 
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This is the function used for the 



fits to the experimental data in Sec. Ill 



III. EXPERIMENTAL & RESULTS 

The diffraction experiments were performed on a Si02 
nanoparticle-polymer grating [36] with a grating spacing 




FIG. 4. Dispersion surface in the off-Bragg geometry. Top: 
overview including the four permitted wave vectors q^s^ ■ 
The orange lines represent part of the dispersion surface. An 
approximate version of the latter - as usually assumed in liter- 
ature - is shown in dotted gray line style. Bottom: A closeup 
of the dispersion surface in the vicinity of the Bragg condi- 
tion with the mismatch Aq. Aqo is the difference between 
the gray asymptotic lines (viz spheres) along the iV-direction, 
i.e., horizontal. 



A = 2tt/\K\ = 500 nm using an s-polarized He-Ne laser 
beam at a wavelength A of 632.8 nm. The sample is the 
same as used for the experiments concerning the direc- 
tion of the diffracted beam discussed in Ref. [27]. A 
schematic of the setup is shown in Fig. [5] The beam 
was prepared using a standard beam expansion system 
(microscope objective lens, pinhole, magnification lens) 
followed by a diaphragm with a diameter of about 5 
mm. The sample is placed on a rotation stage (resolu- 
tion 1/1000 deg) and is rotated around an axis (vertical) 
parallel to the y— axis while measuring the zero and first 
order diffracted powers, i.e., rocking curves are recorded. 
To get reliable values in the far off-Bragg regime we mea- 
sured the diffraction powers using Si-photodiodes which 
were plugged into a light amplifier allowing linear am- 
plification in the 10 5 -order signal range. To ensure most 
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FIG. 5. Experimental setup scheme: L,P,BE,D,BS,S,PD de- 
note the He-Ne laser, polarizer, beam expander, diaphragm, 
beamsplitter, sample and photodiodes, respectively. The am- 
plifier is symbolized by the triangle (top). Sample geometry 
for measuring the diffraction efficiency (left bottom) and the 
background (right bottom), rotation axis was always y. 



reliable results we took at least five power values at each 
angular position thus being able to evaluate a standard 
error of the mean. Furthermore, a background measure- 
ment was conducted by first rotating the sample by 90° 
around the sample surface normal, and then recording a 
rocking curve again, this time the grating vector K being 
parallel to the rotation axis y. By attaching an opaque 
circular mask (diameter 5 mm) to the sample's front- 
surface, it was assured that the incident beam passed the 
identical sample volume for both, the diffraction experi- 
ment and the background measurement. The results of 
the diffraction experiments are shown in Fig. [6] A simul- 
taneous least-squares fit to the weighted and background- 
corrected diffraction efficiency for the minus and plus first 
order were performed in the angular range \0 — 0b\ ^ 0.15 
radian, therefore including the far off-Bragg region with 
18 side minima. It is obvious from the plot that while 
the KVCM does not fit the experimental data, BVM and 
DDT do. This is quantified by the Xk vame being nearly 
an order of magnitude larger than other approaches. The 
most striking difference is observed in the positions of the 
minima. In Fig. [7] the experimental minimum positions 
are compared to those obtained from the theories. From 
the angular dependence of the residuals we conclude that 
the KVCM strongly deviates in the far off-Bragg region. 



IV. DISCUSSION 

As already said in the introduction, typically the dis- 
persion surface in the DDT is approximated in the vicin- 



ity of the Bragg condition so that the dispersion surfaces 
are hyperbolae. This of course is not applicable for the 
far off-Bragg region. On the other hand, a direct com- 
parison of analytical formulae is only possible if proper 
approximations are applied, which is done in Appendix 

El 

Let us start our discussion with mentioning that Aq 
(appearing in DDT) is equivalent to Aqo (appearing in 
the BVM) in the limit of zero coupling (k = 0). In the 
wavevector diagram (Fig. [3| a background refractive in- 
dex no is implicitely assumed to describe the permitted 
wavevectors of propagation in the medium. On the other 
hand, we know that at the exact Bragg condition a wave 
experiences a different refractive index within the limits 
no ± Tii. This is only reflected in the DDT, where - in 
the frame of the two mode case - two waves with slightly 
different refractive indices propagate towards each of the 
two reciprocal lattice points and K shown in Fig. [4] 
The value of Ag is the difference between the two cross- 
ing asymptotic gray lines (i.e. actually circles) along the 
sample surface normal shown in the closeup of Fig. [4j 

In Ref. [27] the effective thickness of the grating - for 
the same sample - was estimated to be do = 58.5±0.05/im 
by fitting Eq. (10) to the angular dependence of the 
diffraction efficiency shown in Fig. [6j The discrepancy 
to the value determined in the present work originates 
from three facts: (i) here we used experimentally mea- 
sured values weighted by their standard error of the mean 
(no weighting in [27]), (ii) the angular range was ex- 
tended to the far off-Bragg region which has particular 
influence on the obtained thickness via the phase func- 
tion through its minima, and (iii) the background was 
thoroughly subtracted. Browsing through the literature 
of evaluating the refractive-index modulation of volume 
holographic gratings, we find mainly two approaches: a) 
measuring a rocking curve and calculating rt\ and d by 
using the KVCM in the linearized version and b) just 
assuming the grating thickness to be the measured me- 
chanical sample thickness and calculating rt\ from the 
diffraction efficiency at the (supposed to be) Bragg an- 
gle. While it is evident that the latter can only serve as 
a rough estimate, the results of the former are usually 
taken as serious parameters. When fitting of the data 
shown in Fig. [6] (without weighting) is performed on this 
basis, i.e., less thorough, values around the Bragg peak 
contribute the most to the result. This paradoxically 
leads to fitting parameters which are much closer to the 
ones obtained by BVM and DDT (deviations of about 
1% for ni and d). The latter, however, is only due to the 
fact that we have an almost ideal grating very accurately 
described by Eq. ( refeq:grating). 

After all we are left with the question: are we able to 
decide if any of the theories is superior? Our experiment 
and the sample under investigation provide almost ideal 
conditions, i.e., the assumptions about the grating taken 
in the models are fulfilled. We identify the KVCM as less 
reliable due to the significantly higher Xk ^ Xbd ( see 
Fig. [6|. A decision on the basis of x 2 alone to recom- 
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FIG. 6. Dependence of the first order diffraction efficiency on the Off-Bragg parameter <& (see Eq. (|3J> ) for a nanoparticle- 
polymer grating. The plus and minus first order data (symbols) were fitted using either the KVCM (green line), BVM (red 
line) or the DDT (black line). Fitting results including the weighted % 2 are given in the inset. For further details, see text. 
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FIG. 7. Comparison of the minimum positions between the 
experimental data and the three models. The deviation of 
the measured minimum positions from the three fitted models 
(residuals) is shown. 



mend either the BVM or the DDT is not evident. How- 
ever, when it comes to less perfect gratings, e.g., taper 
with an attenuation of the grating modulation along the 
sample depth [37] due, for example, to non-neglible linear 



absorption during recording [38], there is no known way 
to apply the modal approach that leads to the DDT. In 
contrast, it has been successfully proven that the prob- 
lem can even be analytically treated using a coupled wave 
analysis [T2] , 



V. CONCLUSION 



Diffraction by a one dimensional volume holographic 
grating with a relatively thick film (about 57 micron) and 
high refractive- index modulation (ni about 0.003) allows 
to demonstrate a significant difference between diffrac- 
tion theories. For now we conclude, that both coupled 
wave approaches, KVCM and BVM, and the correspond- 
ing differential equations are approximations to the DDT. 
While the KVCM fails to give the correct phase function 
for the diffraction efficiency, the BVM differs by a factor 
cr/cs in the amplitude. We finally comment that from 
a practical point of view most of the gratings are not 
that perfect, and thus the BVM including the attenua- 
tion along sample depth as discussed in Ref. [12] should 
be the first choice. 
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Appendix A: Derivation of the dispersion surface 
Eq. (ITil) 



where K g = gK,g £ N. Further, it is known that so- 
lutions to the wave equation for periodic media have a 
particular form 



m^ikx 



(A2) 



We follow the derivations given by Batterman and Cole 
for x-rays [19] and Russell for light [22] . 

We start with the Helmholtz equation (J2|. The refrac- 
tive index from Eq. is a real quantity and periodic 
with the reciprocal lattice vector K = (K, 0, 0); the rela- 
tive permittivity of the medium is defined as usual, i.e., 
n 2 = e e R. Hence, 

e(x)=so + \ + (Al) 

geN,g>0 I 



where S^ 1 are functions with the periodicity of the 
medium. The E m (x) are eigen-modes of the wave equa- 
tion labelled by m (indicating different energies num- 
bered in increasing ord er) an d ar e called Bloch waves. 
By inserting Eqs. ( |AT| and flA2| into Eq. Q we end 
up with an infinite number of algebraic equations for the 
infinite number of coefficients S^P". 



E 



geN,g>0 



E 



(A3) 



This is effectively the Fourier transformed version of the 
wave equation in a periodic medium. 

Assuming that only two of the coefficients considerably 
differ from zero - say - the system of equations 

reduces to 



(P 2 -\tf)S? + ^8,8^ = for k = q 



^0 _ cm 



3 2 - \q + K\ 2 )£™^ = for k = q + K 



Nontrivial solutions to this system of equations require 
the determinant of the coefficient-matrix i.e., the 

k 

characteristic polynomial, to vanish. This leads to a 
quartic equation in \q\ 



W 2 -\ q r)W 2 -W+K\ 2 )=[fe 1 



(A4) 



which gives the moduli of the m permitted wavevec- 
tors for propagation in a periodic medium. The latter 
is Eq. (14) when using the relations n = fco^i/2,£i = 
2nino and the notation q = <f#, q + K — qs . 

Finally, we find the total field in the grating as a sum 



over all eigen-modes 



£%e lk *, (A5) 



where the coefficients v m must be determined considering 
the boundary conditions. We express the ratio 



2k 



cm 



(A6) 



q+K 

and assume that the incident wave at the entrance bound- 



ary has an amplitude equal to unity so that Eq. (A5) 
yields conditional equations for the v m : 



8 1 - 

q+K 



S 1 .£2-£i£° 

q+K q q q+K 



8° - 

q+K 

vi = - v 0gi — • 



q+K 

The field- amplitude of the waves travelling towards the 
reciprocal lattice point K is 

= voS^ge^+XV + v^ge^+X)* (A7) 
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so that for the geometry used in the experiment dis- 
cussed here (Fig. [5| the diffracted amplitude E-i(z = d) 
amounts to 



E-i 



(|g (0) | 2 -^ 2 )(|g (1) | 2 -/? 2 ) 
2/3n(\qW\ 2 - \q(°)\ 2 ) 



(A8) 



Using the definition r\ = \E-\{z = d)/l\ 2 we end up with 
Eq. pi. 



Appendix B: 'Proper' approximation in the DDT for 
comparison with KVCM and BVM 

To obtain an equation similar to the form of Eqs. (|6| 
and (10) also for the DDT, we approximate Aq using Eq. 
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V^ 2 + (2k) 2 
2(tf + /3cos 2 0) 
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Then the functional form of the diffraction efficiency - 
an amplitude term A multiplied by an oscillatory phase 
term sin 2 <1> - looks similar to that of Eqs. ^ and (10) 
yielding 
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(B3) 



(B4) 



For the experimental angular range ##-0.155. . . Ob +0.166 
radian (angles in the medium) and k = 2.08 x 10 _2 /im _1 
(as obtained from the fit) the deviation of Eq. ( |B2|)j rom 
Eq. ( [T6| ) is less than about 0.3%. Rewriting Eq. (|6| for 
comparison in terms of X, vk^k, we have 



A K 



A K sin 2 § K 



(B5) 
(B6) 
(B7) 



& D y/T+X. 

Thus, the phase term of the KVCM differs from that of 
the DDT by a factor \/l + X while the amplitudes are 
identical. We proceed with Eq. (10) in terms of X which 
yields 

T] B = A B sm 2 $ B (B8) 

Ab = 7TTW 
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FIG . 8. for each approach according to Eqs. 

( |B4[ ),(B7),( B10| ). Roots of the diffraction efficiency occur at 
angles for which <3>(0 S ) = sty, i.e, where the dotted horizontal 
lines intersect the functions. 




(2k) 



(/?4(i 



VT+2X)) 2 VT+2X 



(B9) 



VT+2X)) 2 + 



(2k) 



VT+2X 



(B10) 



Despite the formal similarity of the expressions Eqs. (|6| 
and (10) for KVCM and BVM, respectively, we need to 



take a further approximation for a direct comparison. We 
recall that we are interested in the far off-Bragg region, 
where the theories might significantly diff er on e anoth er. 
Then i? > (2k) and we neglect k in Eqs. (TB2} - QBlpp in 



comparison to $ and X. The approximate phase func- 
tions then read 
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\X\ 



1 fic R d 
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II - V1 + 2X 



0,K 
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\x\ 
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(Bll) 
(B12) 
(B13) 



Higher order minimum positions are given by <&(0 S ) = sir 
with s G N. Fig. 8] shows Q>(0) for each of the theories. 
BVM and DDT lead to identical results in the phase 
function up to second order in — whereas KVCM is 
equivalent only in the very vicinity of the Bragg condi- 
tion. For the amplitude part of the diffraction efficiency, 
the situation is different. Assuming again d ^> 2k (for 
the far off-Bragg angular range) we get 

Ao.D = ! 77^7 ) (HID 



A ,B 



(3c 2 R X 
1 



(2«) 



■ 2X (04 [1 



V1 + 2X])V1 + 2X 
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FIG. 9. The amplitu de fu nctions A(6) - Eqs. (|B3| , (|B6J , (|B9 ) 
- and Ao(0) - Eqs. (B14) - in the off-Bragg region for each 
model. 



In Fig. [9] the amplitude functions A(8) are shown. This 
time the DDT is identical to the KVCM and differs from 
BVM, namely by a factor y/l + 2X. The relative devia- 
tion of 1 — ^/T+~2XAb/Ad is less than 0.7% in the range 
discussed. 

To summarize: For the diffraction efficiency 



• the phase term of the DDT agrees to that of the 
BVM excellently; the KVCM differs by a factor of 
(l + X)- 1 ^ while 



• the amplitude term of the DDT is equal to that of 
KVCM; whereas now that of the BVM differs by a 
factor of (l + 2X) 1 /2. 



While the first statement is fully confirmed by the ex- 
perimental data and can be seen in Figs. [6] and [7| the 
second statement can hardly be verified from the dataset. 
One reason could be that the side wings of the amplitude 
function are also influenced by absorption, an issue not 
considered here. 



Approximate theories and their interrelation 

As said above, the equation for the dispersion sur- 
face Eq. (14) discussed in nearly any publication is ap- 



proximated to form hyperbolic sheets in the vicinity of 
the Bragg incidence, i.e., the quartic equation in f3 is 
transformed in a quadratic one with the approximation 
f3 2 - (q) 2 w {P - \q\)2f3. Then the magnitudes of the 
permitted wavevectors are given by 



Am) 



(B15) 



The approximate version of the dispersion surface is 
shown in Fig. [4] as dotted grey lines. Such an approxi- 
mation leads to the identical diffraction efficiency for the 
DDT and the KVCM as given by Eq. Q [22}[39]. Using 
the approximation of Eq. (Bl ) and truncating the expan- 
sion of the p hase correction factor (1+X) -1 / 2 « 1 + 0[X] 



in Eq. (B2 ) with the constant term, we arrive at the same 
result. 

In addition, the off-Bragg parameter is frequently lin- 
earized as in Ref. [TT], so that #(0) = K cos 6b{0 — b) 
which then should give even worse results and anyhow is 
valid only close to the Bragg angle b- 
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